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Abstract 

We determine the critical value of the coupling where the first order quantum phase 
transition takes place for lattice SU( 2) Yang-Mills theories in dimensions higher than four. 
Within a Mean-Field approach we derive an approximate law valid for any dimension d 
and in the context of a Monte Carlo approach, in addition to the already known d — 5 
case, we look at d — 6, 7, 8. 
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Even though high (d > 4) dimensional Yang-Mills theories are perturbatively non- 
renormalizable, one can not exclude the possibility that there exists a regime in their 
phase diagram where a physically useful cut-off effective theory can be constructed. Also 
from the theoretical point of view, these theories contain the basic ingredients (the gauge 
fields) of various held theories motivated by string theory. 1 Therefore, their study is 
potentially useful. 

A possible universal property of high dimensional Yang-Mills theories that emerges 
from previous studies comes from the fact that in five dimensions a ’’bulk” or ’’quantum” 
phase transition appears dividing the confined phase from a Coulomb phase [2], Recall 
that 4 d pure Yang-Mills (for SU(N) with N < 4) 2 at zero temperature has only a confined 
phase where the gauge fields always form flux tubes. Intuition says that as the number of 
dimensions transverse to the surface of the tube increases, the harder is to sustain a stable 
flux tube thus requiring a stronger coupling. In phase diagram terminology, we expect 
the confined phase to persist but also to shrink as d increases. Here we will perform a 
first check of this statement. 

We take SU (2) as our model and regularize it on a d-dimensional, Euclidean, periodic 
lattice with lattice spacing a and linear dimension L. The action is the standard Wilson 
plaquette action 

s e = P'%2 l -\ il i U A x )} »■ (!) 

x l<Mi v<d L 

where U liu (x) = U ll (x)U u (x + ajijUjJx + au)Ul(x), /q v — 1, • • • , d is the elementary pla¬ 
quette located at the site x with U^{x) = e mA U x ) and 6 — A) represents the dimensionless 
lattice coupling. We will employ two methods, on one hand the Mean-Field (MF) approx¬ 
imation [4, 5], an analytic method expected to work well near the phase transition and in 
general increasingly well as d grows 3 and Monte Carlo (MC) simulations on the other. 

For SU{ 2) in d dimensions, the (ungauge-fixed) MF approach to zeroth order deter¬ 
mines the confined and Coulomb phases via the solution to the coupled equations for the 
MF background v 0 [4] 

Vo = ho = 2vjftd - 1 )j3 (2) 

J i(no ) 

with I„(ho) the modified Bessel function, by defining the Coulomb phase as the regime 
of (3 where there is a solution with no ^ 0 and as the confined phase otherwise. In [5] 

the equations were solved for d = 5 by a numerical, iterative method. The smallest 

positive and real non-vanishing value of the background vq c satisfying eqs.(2) (i.e. where 

1 For Monte Carlo studies of high dimensional (supersymmetric) Yang-Mills theories from the point of 
view of matrix models, see [1]. 

2 For TV > 4 a first order bulk phase transition emerges already in d = 4 [3]. 

3 Or when N of SU(N) grows; this however tends to shrink the Coulomb phase instead [6]! 
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the iteration stabilizes), determines the critical value of the lattice coupling f3 c where the 
phase transition takes place. It is expected to be a quantum phase transition since the 
MF at this order is volume independent. This of course needs to be checked. In fact, 
for d — 5 it was found by a MC simulation on a 4 5 lattice in 1979 by Creutz to be a 
quantum, first order phase transition. Subsequently this was confirmed (and extended to 
anisotropic lattices) by several authors [7]. Apart from the fact that both methods agree 
on the order of the transition, their estimates for the value of the critical coupling are also 
quite close: the /3^ F ~ 1.6762017 of the MF (corresponding to v 0c — 0.73333) [5] is to be 
compared with the (3^ c = 1.642 ± 0.015 of the MC [2], 

An observation stemming from eqs.(2) is that the quantity B = (d — 1)/^ F and 
therefore also vq c are d- independent. We can then solve eqs.(2) for general d by noting 
that the zero of the function F = / 2 (/io)/ Ii(ho) — v 0 that signals the phase transition is 
one where F(v o) has an extremum. The extremization condition F' — 0 yields 


3Bi>l 



1 . 


(3) 


Using eq. (2) and the recursion identity I v {ho) = I u - 2 (ho) — 2(y h ^ J 1/ _i(/i 0 ) we can express 
I3/I1 = 1 — 4u 0 /ho and Io/h — 2//t 0 + Vo and reduce eq. (3) to the quadratic equation 

x 2 ~ x + ^ = 0 (4) 

with x = Vq. The solution determines the value of the background at the critical point 


1 // 20 1 

= 71V v 1 " Tb ' (5) 

Substituting the above back in the equation F = 0 results in an algebraic expression with 
only parameter B , whose relevant root can be found numerically to be B ~ 6.704840, 
determining v 0c ~ 0.7333 from the upper sign of eq.(5), as expected. Thus, we find that 
the equation 

(d - 1)/? C MF ~ 6.704840 (6) 

fixes the SU (2) critical coupling in any dimension d > 4. 


What makes it possible to go high in d with Monte Carlo simulations is that we are 
dealing with a bulk phase transition. This means that as long as the lattice extent is 
large enough so that finite size effects do not interfere, the phase transition is visible. 
Most times a A d lattice will suffice to observe the effect, even though larger lattices will 
be clearly needed to describe it with better precision. The order parameter used in order 
to locate the phase transition is the plaquette 


P 


d(d- 1)L D 


£ £ 


x 1 <fj,<u<d L 


-tr {U^ v {x)} 


(7) 
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The Kennedy-Pendleton heat bath algorithm [8] has been used to update the gauge field. 
Far from the region of the phase transition the number of thermalization sweeps needed 
so that the plaquette achieves its equilibrium value has been of the order of 100, 150, 200 
for d — 6, 7, 8 respectively. We have also employed overrelaxation hits to decorrelate the 
measurements. They don’t seem to affect the results very much, since the measurements 
of the autocorrelation times yield very small values. The acceptance rates have been of 
the order of 70 — 90 per cent. 

The phase diagram has been obtained for d = 6, 7, 8, while we also quote the results 
for d — 5. The lattice sizes used have a linear dimension L — 4 and, after thermalization, 
hysteresis loops have been performed. The step in fd was 0.07 (starting at fd = 0.40 and 
going up to fd = 1.8 and back); after 1000 initial heat up sweeps, 8000 iterations through 
the lattice have been done at each fd value. There exist well known approximations that 
we have used as guides. In the strong coupling regime (small values of fd) the plaquette 
P is well approximated by P ~ |, while in the weak coupling the approximations reads: 
P ~ 1 — Jk. It appears that the fd values shown in the figure are not large enough 
to approach the weak coupling limit, but we have checked that for larger values the 
agreement is good. The results for the hysteresis runs along with the strong and weak 
coupling approximations are shown on the left in Fig. 1. The phase transition is seen to 
be a strong first order one, even for this quite modest linear dimension of the lattices. We 
can read off estimates for the [d intervals in which the critical values lie: 1.29 < fd < 1.45 
for d = 6, 1.07 < fd < 1.27 for d = 7, and 0.86 < fd < 1.04 for d = 8. The right of Fig. 1 
depicts a comparison of the mean field estimates for the critical fd and the corresponding 
analytical expression (6) against the results of the Monte Carlo runs. We have also quoted 
the five-dimensional result from [2] for completeness. We observe a quite good agreement 
between the two methods. 

The determination of the order of the phase transition as well as of the pseudocritical 
values for fd is usually done through long runs and by examining the fluctuations of 
the plaquette between the values pertaining to each of the two metastable states. This 
permits the construction of histograms and the determination of the volume dependence 
of the specific heats. Unfortunately, the gap between the two values grows so much with 
increasing dimensions, that it has proven impossible to observe this behavior with the 
algorithm used. This would require a multicanonical simulation, so we defer this to a 
forthcoming publication. For the time being, we have used long runs (of the order of 10 5 
iterations) just to be sure about the upper and lower bounds that we report. The results 
are encoded in the upper and lower bounds for the critical fd reported above. 

Extending the work of [2] we determined the critical value of the coupling where a first 
order bulk phase transition takes place for high dimensional SU (2) lattice gauge theories. 
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Figure 1: Left: Hysteresis loops for 6, 7 and 8 dimensions. The strong and weak coupling 
predictions are also included. Right: The pseudocritical values for /3^ c versus and their 
estimated errors. The Mean-Field prediction eq.(6) is represented by the straight line. 


We first derived a law valid in any dimension d > 4 based on the Mean-Field method and 
then performed corresponding Monte Carlo checks for the first time in d = 6, 7 and 8 
dimensions. 

The values of the plaquettes at selected values of /3, as well as several results relating 
to the error and autocorrelation issues are presented in the Appendix. 
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1 Appendix 

We collect in this Appendix some details of the simulations. Near the region of the 
phase transition the autocorrelation question becomes important. For the six-dimensional 
models we plot in Figs. 2 the mean plaquette resulting in long runs for a hot and a cold 
start and values of (3 smaller than the phase transition value. We see that the hot start 
makes the plaquettes take on their final values very soon, while the cold start spends 
considerable computer time in a false vacuum before actually landing on its true value. 
The time spent in the false vacuum becomes longer as one approaches the phase transition. 
For the relatively large (3's depicted in Figs. 3 the opposite effect takes place: it is the hot 
start which gives a fluctuation around the false vacuum and this fact is more intense as one 
approaches the phase transition region; on the contrary, the cold start fluctuates around 
the true vacuum from the beginning. It appears that the measurements in this region 
are highly correlated; however the correlation functions indicate that the autocorrelation 
time is bigger than the total computer time of the run, so it makes no sense to display 
them. This behaviour just signals the limitations of the method. Similar results for d = 7 
are plotted in Figs. 4 and 5. Similar behaviour is observed for d = 8, but one has to 
fine tune too much to achieve the corresponding wandering around the false vacua, so we 
don’t show the results. 

It is useful to report more details on the values of the plaquettes and their errors, 
as calculated by the Jackknife method. We report the relevant values at the values of 
f3 that have been used for the determination of the phase diagram in tables 1, 2, 3 for 
d — 6, d — 7, d — 8 respectively. We also depict the results of the long runs in Fig. 
6, including the error bars. We show both hot starts and cold starts along with the 
corresponding errors. The errors grow large only near the phase transition, as expected. 
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Time evolution 

P = 1.289 P'g = 1.289 iterations = 100000 


Time evolution 

P g = 1.290 P'g = 1.290 iterations = 100000 


Time evolution 

P = 1.291 P' = 1.291 iterations = 100000 
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Figure 2: Long runs for six dimensions below the phase transition, at fi = 1.289, 1.290 and 1.291. 
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Time evolution 

P = 1.427 p' = 1.427 iterations = 100000 


Time evolution 

Pg = 1.426 p'g = 1.426 iterations = 100000 


Time evolution 

pg = 1.425 P'g = 1.425 iterations = 100000 
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Figure 3: Long runs for six dimensions above the phase transition, at P = 1.427, 1.426 and 1.425 







































Time evolution 

Pg = 1.0728 P'g = 1.0728 iterations = 100000 


Time evolution 

P g = 1.073 P'g = 1.073 iterations = 100000 


Time evolution 

Pg = 1.0732 P'g = 1.0732 iterations = 100000 




Figure 4: Long runs for seven dimensions below the phase transition, at (3 


1.0728, 1.0730 and 1.0732. 


Time evolution 

Pg = 1.3105 P'g = 1.3105 iterations = 100000 


Time evolution 

Pg = 1.310 P'g = 1.310 iterations = 100000 


Time evolution 

Pg = 1.3095 P'g = 1.3095 iterations = 100000 
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Figure 5: Long runs for seven dimensions above the phase transition, at (3 = 1.3105, 1.3100 and 1.3095. 
































Figure 6: Plaquette values with errors in 6,7 and 8 dimensions 
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Table 1: Plaquette values & errors in 6 dimensions 


b 9 

Cold start 

Cold start 

Hot start 

Hot start 

0.8500 

0.210118 

0.000003 

0.210116 

0.000005 

0.8900 

0.220185 

0.000005 

0.220185 

0.000003 

1.0700 

0.267476 

0.000003 

0.267473 

0.000005 

1.0728 

0.275 

0.007 

0.268251 

0.000003 

1.0730 

0.30 

0.02 

0.268303 

0.000006 

1.0731 

0.290 

0.018 

0.268336 

0.000003 

1.0732 

0.30 

0.02 

0.268388 

0.000009 

1.0733 

0.285 

0.017 

0.268393 

0.000006 

1.0734 

0.30 

0.02 

0.268415 

0.000003 

1.0735 

0.437191 

0.000017 

0.268445 

0.000003 

1.1680 

0.540070 

0.000003 

0.296037 

0.000004 

1.1720 

0.542651 

0.000004 

0.297304 

0.000004 

1.1750 

0.544582 

0.000003 

0.298252 

0.000004 

1.1980 

0.558431 

0.000005 

0.305641 

0.000005 

1.2800 

0.599198 

0.000005 

0.336737 

0.000008 

1.3090 

0.611366 

0.000004 

0.59 

0.03 

1.3095 

0.611544 

0.000002 

0.55 

0.04 

1.3100 

0.611763 

0.000002 

0.597 

0.015 

1.3105 

0.611938 

0.000003 

0.609 

0.002 

1.3120 

0.612546 

0.000004 

0.612549 

0.000004 

1.3300 

0.619573 

0.000003 

0.619576 

0.000004 

1.3400 

0.623345 

0.000004 

0.623346 

0.000004 


b 9 

Cold start 

values 

Cold start 

errors 

Hot start 

values 

Hot start 

errors 

0.8400 

0.208402 

0.000001 

0.208402 

0.000001 

0.8900 

0.221268 

0.000001 

0.221268 

0.000001 

0.9000 

0.223868 

0.000003 

0.223873 

0.000003 

0.9100 

0.226487 

0.000001 

0.226487 

0.000001 

0.9200 

0.229113 

0.000004 

0.229110 

0.000002 

0.9215 

0.2301 

0.0006 

0.229511 

0.000004 

0.9350 

0.451421 

0.000010 

0.233080 

0.000003 

0.9400 

0.459860 

0.000006 

0.234412 

0.000002 

0.9450 

0.467337 

0.000011 

0.235743 

0.000004 

1.0000 

0.524477 

0.000005 

0.250721 

0.000003 

1.0500 

0.560293 

0.000004 

0.264967 

0.000003 

1.1000 

0.589056 

0.000004 

0.280106 

0.000001 

1.2000 

0.634465 

0.000001 

0.316682 

0.000002 

1.2100 

0.638365 

0.000006 

0.321559 

0.000005 

1.3000 

0.669492 

0.000003 

0.669492 

0.000003 
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0.697810 

0.000004 

0.697810 

0.000004 


Table 3: Plaquette values & errors in 8 dimensions 


Table 2: Plaquette values & errors in 7 dimensions 























































































